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We calculate the coupling constants between the light vector mesons and heavy mesons within the 
framework of the light-cone QCD sum rule in the leading order of heavy quark effective theory. The 
sum rules are very stable with the variations of the Borel parameter and the continuum threshold. 
The extracted couplings will be useful in the study of the possible heavy meson molecular states. 
They may also helpful in the interpretation of the proximity of X(3872), Y(4260) and Z(4430) to 
the threshold of two charmed mesons through the couple-channel mechanism. 
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I. INTRODUCTION 



\q ' A number of hadronic states which can not be easily accommodated in the conventional quark model have been 
observed experiementally in r ecent years, such_as A(3872) T], Y(4260) and Z+(4430) 0. Their masses are very 
close to the thresholds of DD* , D*D* and D*D 1 respectively. It was speculated that the coupled-channel effect may 
play an important role because of the attraction between these D mesons. Alternatively, they were considered to be 



possible candidates of the heavy molecular states composed of two D mesons. These loosely bound states are formed 
Qh) by exchanging light mesons such as ir, a, p and u> etc. Up to now, the pion heavy meson strong interaction is relatively 
known due to chiral symmetry. However, the vector meson heavy meson strong interaction has not be extensively 
studied yet, which accounts for the relatively short distance interaction between two heavy mesons. 

Heavy quark effective theory (HQET) 4j is a systematic approach to study the spectra and transition amplitudes 
of heavy hadrons. In HQET, the expansion is performed in term of 1/rriQ, where tuq is the mass of the heavy quark 
involved. In the limit itiq — > oo, heavy hadrons form a series of degenerate doublets due to heavy quark symmetry. 
The two states in a doublet share the same quantum number ji, the angular momentum of the light components. The 
B(D) meson doublets (0~, 1 _ ), (0 + , 1 + ) and (1 + , 2+) are conventionally denoted as H, S and T. 

Light-cone QCD sum rules (LCQSR) [5| is a very useful non-perturbative approach to determine various hadronic 
transition form factors. One considers the T-product of the two interpolating currents sandwiched between the vacuum 
and an hadronic state in this framework. Now the operator product expansion (OPE) is performed near the light-cone 
rather than at small distance as in the conventional SVZ sum rules |6j. The double Borel transformation is always 
, invoked to suppress the excited state and the continuum contribution. 

The p coupling constant between D and D* was calculated with LCQSR in full QCD in Ref. |7j. The couplings 
9h*h*p, Jh*h*p, QHHp and fa*Hp were calculated in full QCD in Ref. Q. Their values in the limit tuq — + oo are 
also discussed in this paper. The p coupling between doublets T and H are studied in the leading order of HQET in 
Ref. 1 . 

In this work we use LCQSR to calculate the p coupling constants between three doublets H , 5, T and within the 
two doublets H , S. Due to the covariant derivative in the interpolating currents of T doublet, the contribution from 
the 3-particle light-cone distribution amplitudes of the p meson has to be included when dealing with the p decay 
between doublets T and H(S). We work in HQET to differentiate the two states with the same J p value yet quite 
different decay widths. The interpolating currents J"p adopted in our work have been properly constructed in 
Ref. [nj. They satisfy 

(0\J£j?(P)\3',P'Jl) = fPnSjj'SpP'S^-^, (1) 

1 (o|t{j^(x)j;^'(o)}|o) = s jjf s PP ,6 jlf (-iysgr Pl ---g^ f dts(x-vt)n Ptjl (x), (2) 
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in the limit itiq — > oo. Here r] ai ' aj is the polarization tensor for the spin j state, v is the velocity of the heavy 
quark, g" 13 = g a/3 — v a v@ , S denotes symmetrizing the indices and subtracting the trace terms separately in the sets 
(a i ■ ■ ■ otj) and (/3i • • • (3j). 



II. SUM RULES FOR THE p COUPLING CONSTANTS 

We shall perform the calculation to the leading order of HQET. According to Ref. the interpolating currents 
for doublets H, S and T read as 



Jl 



1,+, 



J, 



tc*ia2 
2,+,l 



h v q, 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



where /i„ is the heavy quark field in HQET, 7^ = 7^ — th)**, 2?^ = — (T> ■ v)v^ : = g^ v — v^v v ', and is the 
velocity of the heavy quark. 

We consider the p decay of T\ to H\ to illustrate our calculation. Here the subscript of T(H) indicates the spin of 
the meson involved. Owing to the conservation of the angular momentum of light components in the limit itiq — > 00, 
there are three independent p coupling constants between doublets T and H . We denote them as g^ Hp , g'rHp an< ^ 
d^THp where s,p,d - ■ ■ and the number following them indicates the orbital and total angular momentum (I, jh) of the 
final p meson respectively. All of these three coupling constants appear in the decay process under consideration. The 
decay amplitude can now be written as 



n dl 



n d2 



(9) 



where 77, e* and e* denote the polarization vector of Ti, Hi and p respectively, q is the momentum of the p meson, 
q 2 = m 2 and gf = q^ — (q ■ v)v^. 1=1, l/\/2 for the charged and neutral p meson respectively. The vector notations 
in Levi-Civita tensor come from an index contraction between Levi-Civita tensor and the vectors, for example, 

To obtain the sum rules for the coupling constants gf\# lP , QtxHxp anl ^ 9t\h 1 o' we consider the correlation functions 



d 4 xe- ik - x (p(q)\T{j[_ h (0)Jl^(x)}\0) = Ii\e^G^ lHlp (uj,u;') 

e af3qv {e* -q t )-\e al3e ' v q 2 t 



e *e* q v ql 3 + e /3e*qv q a 



where w = 2v ■ k, u>' = 2v ■ (k — q). In the leading order of HQET, the heavy quark propagator reads as 

1+v 



(0\T{h v (0)h v (x)}\0) 



dt6 4 (-x~vt). 



(10) 



(11) 
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The correlation function can now be expressed as 



^ I dxe- ik * I™ dtS(-x - vt)Tv j 7 f i±H(-i 75 )(D« - l>y*T> t )(p(q)\ q (x)q(0)\0) }■ <l->> 



It can be further calculated using the light cone wave functions of the p meson. To our approximation, we need the 
two and three-particle light-cone wave functions. Their definitions are collected in the Appendix [Bl 
At the hadron level, the G's in (fT0|) has the following pole terms 



G Ti h iP (uj,uj ) = —= —= - + -= - + -- , (13) 

(2A_ )1 /2-w')(2A +j3 /2-a;) 2A_ il/2 - w' 2A +i3/2 ~uj 

where A 1/2 = toh — mg, A + 3 / 2 = — toq. f-,1/2 etc is the overlap amplitudes of their interpolating currents 
with the heavy mesons. 

GtiHip('j- 1 , u>') can now be expressed by the p meson light-cone wave functions. After the Wick rotation and the 
double Borel transformation with w and a/, the single-pole terms in (| 13[) are eliminated. We arrive at 

t— , A +,3/2+ A -,l/2 

-±- f T m jA T (u )u ± + \fJm*Bf{u )± - \fim*C™{fr)Tf Q (±) - J/ p T m2^(u )T/ (^) 
-1/Jm>f )(u )T/ (^) + l/J^(^)'(« )T/o(^) + ^m^^T/o^) 
-I/J m ^ )T/ (^) + ±fJ m fr ± (u )u Tf {^) - l/J^4 2 )(^ )T/ (^) 
+ l/Jm^M T )( 2 )(u )r/ (^) - i/Jm^4 1 ](u )T/o(^) + l/j4% )T 3 / 2 (^) 
-^/ p T (^-) (2) («o)T 3 /2(f ) + Ifimpl-Wto)! - pj m ^° ] (u ) I 

+ l/Jm 2 r[^]( Uo )r /o (^) + I/J TO ^°l( W o)T/o(f ) - |^m^-^(«o)i 
-i/>2r 3 ^ ]( M0 )r/ (^) + -L/^ A[ 2 ]( - o) ^_ _ _L /pm ^[i] ( - o) ^_ + }_ fpm s (uA) m {Uo) A_ 

+ \f P rnlcM{u Q )^ - lf p m 5 p CW(u )± + lf p m%uCf\u )^ + ±f p m 3 p A(u ) 
+ lf p m 3 p A'(u ) - L fprn ? p{u A)'{u Q ) - lf P m 3 p A^(u ) + -/^C^o) + -/pm^fa) 
-lf p m 3 p (uC)W{u ) + l/ pTO 3 3i Q)[11 ("o) - ^f P m 3 p 9^(u ) + (u )u 
+\f P m 3 p g^ m (no) - l/^P.°]( Uo )T 2 /i(f ) + ~f p m 3 p V™(u ) + i/^V^M^/iCf ) 
-\f P m 3 p ipf{u Q ) + \f P ml^(uo) - l -f p ml{u n ) [1 \u Q ) + l/^M^-L + l -f p mlV^\u ) 
-\f pm ^\u ) + l /pTO ^[^o] M i_ _ I/ pm 8$M (tto) + i/, ra /)W( So )T 2 /4) 
+ ^/p-P-9i a)(2) (^)T 2 /i(^) - ^ W^M^MT 2 /^) + i/pmpffW^/iC^) 
-^./p«W||(«o)T 2 ACy) - ^/p"1p^|(uo)T 2 /i(y) + ^/ p 77i p (u^||)'(i2o)r 2 /i(^), (14) 

/— Ji A +,3/2+ A -,l/2 

V^U,/-,*/ +)i e ^ 

= -/J m 2^W N ) 1 + ^mJ(«ftj[)W(iio)i + 2 / > 2 5[- 1 ^K)i + 1/J TO 2 A T (S„)1 



l/Jm 2 A T (u )« l - 2/Jm 2 4 3l (w )i - |/pV(So)T/ (^) + I/J^(iz ) Wo (^ 
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+I/ pm 3AW(uo)p - i/ P m»(uA)W(«o)^ - 4/ p m^W(u )^ + 2/ p m 3 C[ 3 ](<2 )^ 2 
-2f p m 3 p (uC)®(u )± - \f p m p A^\u Q ) - ~f p m p g^ m (u Q ) - lf P m p g^(u ) 
+ ^fp m p9±\ ;n -())uo - fpm p g { l )[2 \u Q ) + f p m p V [0 ' 0] (u ) + f p m p (p^ ] (u ) - f p m p Lp^ ] (u ) 
+f p m p (u n p(u ) - 4/,m p * [ - 2 ' 0l M^ - 4/ p m p $[- 2 '°](u )^, (15) 

/ — ,„ A +,3/2 + A -,l/2 

V3g^ lP /_,i/ +je t 

= ^/J m ^T(uo)^ - ^/J^ p ^t(uo)wo^ - ^/Jpi.(wo)77o(^) + |/J<^±(wo)uor/o(^) 

-3fJm 2 p Tt m M^ - 3fJmlTt lfi] (u )^ + \f P m p A^\up) - \f P m p g { l\u ) + \fpmpgf (u )u , (16) 

where f n {x) — 1 — e x Sfc=o TT ^ s ^ ne con tinuum subtraction factor, and oj c is the continuum threshold, uq — T ^ T , 
T = Ti+t-2 an< ^ = 1 ~ u o- 2\ and T 2 arc the two Borel parameters. We have used the Borel transformation 
B^e auJ = 5 (a - 5;) to obtain ([13 )) . ([T5 J) and (JTBJ). In the above expressions we have used the following functions 
F [a] (u ) and ^[°' 6 ](tio). They are defined as 

pu Q rxz rx 2 

T [n] {up) = / •••/ / T{x 1 )dx 1 dx 2 ---dx n , (17) 
Jo Jo Jo 

^(«o) ee r Z 1 ^ 2 J7(1 ~ a2 " a3 ' a! " a3) d a 3da 2 , (18) 

JO Ju -a 2 a 3 



(20) 



/o uo — a 2 J a 3 



"3 



Ct3=f0— "2 



1 - M0 , d[.F(l-a 2 -a 3 ,a 2 ,a 3 )]/c>a 2 
cta 3 

/0 "3 

^l -1 ' ! (tt ) = / / J"(l - a 2 - a 3 , a 2 , a 3 )da 3 da 2 

+ " p K-. 2m i-a 2 -a 3 , Q2 ,a 3 ) da3dQ2; (21) 

JO Ju -a 2 a 3 

f«0 f»0-t«2 f«3 

J r ' _2 '°'(uo) =11 / .F(l — a 2 — a 2 , x)dxda 3 da 2 

Jo Jo Jo 

2 f«o puo—a 2 

+ - / / «3-? r (l - «2 - CK3, a 2 , a 3 )da 3 da 2 

2 Jo Jo 

1 r° r 1 -* 2 (up - a 2 f 

+- / / .F(l - a 2 - a 3 , a 2 , a 3 )da 3 a!a 2 , (22) 

2 JO Jua-a-2 a 3 
ruo ru -a 2 rot3 rx 2 

.F^ 3 ' ! (tto) =11 / / ^(1 — a 2 — Xi,a 2 ,xi)dxidx2da3da2 

Jo Jo Jo Jo 

^ r u o rUQ—a 2 ftt 3 

+— If I xT{\ — a 2 — x, a 2 , x)dxda 3 da 2 

2 Jo Jo Jo 

Y f u o ruo~a 2 

+ - / / a 3 F(l - a 2 - a 3 , a 2 , a 3 )da 3 da 2 

Jo Jo 

+ - / / .F(l - a 2 - a 3 ,a 2 ,a 3 )da 3 da 2: (23) 

6 Jo Ju - a2 "3 
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III. NUMERICAL ANALYSIS 

In our numerical analysis, we need the mass parameters A's and /'s, the overlapping amplitudes of these interpo- 
lating currents. We adopt A_ 1/2 from Ref. [11]: A_ 1/2 = 0.5 GeV, /_ 1/2 = 0.25 ± 0.04 GeV 3/2 . A +)1/2 , /+,i/ 2 , 
A+,3/2, and /+, 3 / 2 are given in Ref. [lj]: 

A+,i/ 2 = 1.15 GeV, / +)1/2 = -0.40 ± 0.06 GeV 3/2 , 
A +i3/2 = 0.82 GeV, / +j3/2 = 0.19 ± 0.03 GeV 5/2 . 

The parameters appear in the distribution amplitudes of the p meson take the values from Ref. [14]. We use the 
values at the scale fj, = 1 GeV in our calculation under the consideration that the heavy quark behaves almost as a 
spectator of the decay processes in our discussion in the leading order of HQET: 



/ P '[MeV] #[MeV] 


a 2 a 2 Cjp 


&3p w i w 3~p Cl 


4 U 




216(3) 165(9) 


0.15(7) 0.14(6) 0.030(10) 


-0.09(3) 0.15(5) 0.55(25) 0.07(3) 


-0.03(1) -0.03(5) 


-0.08(5) 



We will work at the symmetry point, i.e., T% = T 2 = 2T, uq = 1/2. This comes from the observation that the mass 
differences between H, S and T are less than 0.4 GeV in the leading order of HQET. They are much smaller than 
the Borel parameter T\ , T 2 ~ 3 GeV used below. On the other hand, every reliable sum rule has a working window 
of the Borel parameter T within which the sum rule is insensitive to the variation of T. So it is reasonable to choose 
a common point T\ — T 2 at the overlapping region of T\ and T 2 . Furthermore, choosing T\ = T 2 will enable us to 
subtract the continuum contribution cleanly while the asymmetric choice will lead to the continuum substraction very 
difficult [H[. 

From the requirement that the pole contribution is larger than 60%, we get the upper bound of the Borel parameter. 
This leads to T < 1.7 GeV. The convergence requirement of the operator product expansion leads to the lower bound 
of the Borel parameter T — 1.3 GeV, starting from which the stability of the sum rule develops. The resulting sum 
rules are plotted in FigHH and [3] in the working interval 1.3 GeV < T < 1.7 GeV and uj c = 2.8, 3.0, 3.2 GeV. 

□ 

0.04 



„ 0.03 



< 0.02 
= 
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T[GeV] 

FIG. 1: The sum rule for gT 1 H 1P .f-,i/2.f+,3/2 with u c = 2.8, 3.0, 3.2GeV 

Other p coupling constants between H, S and T doublets can be calculated in the same way as Qt x H\p- Their 
definitions and the relevant correlators are given in Appendix [Al Here we simply present the sum rules for these 
coupling constants: 

JS^Ai - ^e^{- /p m^W(,2 )i,-8/ p m p C[ 3 ]( t 2„)i, 
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FIG. 2: The sum rule for ff$itfip/-,i/2/+,3/2 with <^ c = 2.8, 3.0, 3.2GeV 
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FIG. 3: The sum rule for ffT?H 1P /-,i/2/+,3/2 with u c = 2.8, 3.0, 3.2GeV 



+4/Jm^f 1 (fio)I + *f P m p <p§ ] («„) } , (24) 

< ffip /!,| = ^e^|/Jm^ T (uo)^~2/ p m p .gi a) (u )-4/>x(So)T/o(^)}, (25) 

-8f pm 3 p d 3 \u )^ +4/ /> m^ ] («o)}, (26) 

= ^e^|-/Jm^ T (« )I+4/J^(^)T/o(^)-2/ p m p5 i Q) (,2 )|, (27) 

1 A +|1/2 +A_ |1/2 f 2^[1]/- \ , ,T 2*/ 



-4/,V ± (« )T 2 /i(^) - 8.W^ ) (uo)T/o(^)}, 

-2^ 2 ](Go)/„m^ - 8/ p m p ^ )[2l (S )i + 8/ p m^[ 2] (« )^}, 
v35T 1 s 1(9 /+,i/+,| e 

-l/j m 2 ^(uo) + ^fJm%uA T y(u Q ) - ^/> 2 4 21 M + J/^Wr 2 /^) 

-|/J(^)'(«o)T 2 /i(^) - ^/Jm^M + ^m 2 p Tr ] (uo) A/ pm ^W(« )I 
- /pm ^[-^]( U0 )I + i/^Vl- 1 ' ]^)^ - i/pm^MTM^ ) 
-^/p^p9i a) ("o)T.fo(^) - \f pmp g^ [1 \u )Tf (^) - \f p m p g^ (u )Tf (^) 
+ I/ pTOp3 [% ) Wo (^) + \f P m p V^(u )Tfo(^) + ^f p m p ^ ] (u )Tf (^), 

^ P 1 , , _ A + .3/2+ A +,l/2 

|/Jm 2 4 2l (So) - J/> 2 4%o) + J/Jm 2 ^)! 1 ]^) - ^> 2 ^M 
+ |/> 2 (u^)W(So) + |/> p 4 ] (So)^ - |/> p M T )W(u )^ + ^/Jm^o) 
-4/ p T m 2 M T )'(« ) - |/J m ^W(n )^ + ^/Jm^M^ 
-^/J^B^l^-L - ^/> 2 4 2l ("o) - |/Jm 2 4(So) + j^/Jm 2 ^) N (u ) 

-^/p T ^(-o)t 2 A(^) + ^/J(^)'(So)T 2 A(^) + j|/> P ^ 2 ' 01 M^ - ^/p T ™^ [0 < 0l M 
+y /p T -p^ h2 ' 01 K)^ + 4«rJ M M + f /p T ^ h2 ^ 01 (-o)^ + 4« r i°' 0l M 

-A/ pm 3A(«o)^ + ^/ p m^(S )«o^ - ^/ P m 3 .g^[ 3 ](S )I + jj/ p m^ 2 ](«o)± 
-^/ p m p (^))[ 2 l(fIo)^ + ^/pmpVl-^M^ + |/pm>| 3 l(fIo)^ - \f P ml^{u^ 

+ 6 -f pmp (u n ) l2] (u )± - ^ /pmp ^- 3 -°] M ^ - !/p^ [ - 3 < 01 M^ + |/p-^ h3 - 01 K)^3 

+ ^/p- 3 $ [ - 1 '° 1 K)^ ~ \f P m^\u a )± ^f p rn p9 ( f m (u )TM^) 

+ ^f P m p9 ( f ) (uo)T.U^) ^f P m p9 ^\u )u Tf (^) + lf pmp V^(u )Tf (^) 

+ ^/p m P^I| 11 («o)Tf (^) + ^f p m p tp\\{uQ)T.f Q (^) - ^f p m p y\\{u )uoTf (^), 

ft n C t A +,3/2+ A + ,l/2 

= -3/J^ ] («o) + 3/>^)W(u ) + |/p T ™ 2 AW(So)^ - fjMr) 111 ^)^^ 
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-24/>^4%o)^ + 12/Jm^4 3 l(«o)^ - 12fJml(uB T f\uo)± + e/Jm^-^M^ 

-|/ p m3( u A)[ 2 ](S )^3 - 6/ p m p . 9 ^ [3l (S )l + 6/ p m p ^ )[2l M^ - 6/ p m p (u 3 ^)[ 2 ] («o)± 
-6/ p m p V hl ' 01 (u )^ + 6/ p m^[| ] (!io)i - 6/ p m p <^ 2] (u )^ + 6/ p m p (M^|) [2l (u )^ 
-12/ p m3 V [-3,o ](uo) ^_ _ 24/,m3* [ - 3,01 («o)-L + 24./ p m p *[- 3 >°] («o)-L , (32) 

Due to heavy quark symmetry, the p coupling constants with the same (l,jh) between two doublets are not inde- 
pendent in the leading order of HQET. The values of these coupling constants multiplied by the decay constants of 
the initial and the final heavy mesons are: 

9h h 0P = -~9 P h iHiP = -0-32 ± 0.04 GeV 2 , 

~9h\h 0P = ~9 P h\h iP = -0-46 ± 0.01 GeV 2 , 

~9s h iP = -~9s\h 0P = -9s\ HlP = -0.39 ± 0.03 GeV 3 , 

~9%h iP = -9s lHoP = -~9s lHlP = -0-38 ± 0.06 GeV, 

9%s p = -9&s lP = -OM ±0-03 GeV 2 , 

9 P s\s oP = -9 P s\s ip = -0M±0M GeV2 ' 

3tU„p = - 2 9t[h iP = -S^tU.p = -0.04 ±0.002 GeV 4 , 
~9t iHoP = l~9 d T\ HlP = -2\H~9tIh iP = -0-47 ± 0.06 GeV 2 , 
9t iHiP = \H~9tIh 0P = ^~9t 2 h iP = -0-15 ± 0.01 GeV 2 , 



9t iSoP = 1~9t\s iP = 2y 3#Up = 0.27 ±0.02 GeV , 
~9tIs iP = -\fl~9 P T 2 s 0P = ~^~9 P t 2 s iP = -0-28 ± 0.02 GeV 3 , 

~9 f T lSlP = - Vf ~9 S t\s 0P = -^~9%s lP = 0.31 ± 0.04 GeV, (33) 

where g P H H = 9H H f- 1/2 e ^ c - The errors come from the variations of T and lu c in the working region and the 
central value corresponds to T = 1.5GeV and oj c = 3.0GcV. The g's with their errors are 

9 P h h qP = -9% lHlP = -5-1 ± 0-6 ± 1.3 GeV- 1 , 
&H0P = & xHxP = -7.4 ± 0.2 ± 1.8 GeV" 1 , 
9s\h iP = -9s\ HoP = -9s\ HlP = 3-9 ± 0.3 ± 1.0, 
9% HlP = -9 d s\ HoP = -9 d s\ HlP = 3-8 ± 0.6 ± 0.9 GcV" 2 , 
9 P s° s 0P = -9 P s° lSlP = -2-0 ± 0.2 ± 0.5 GeV" 1 , 
9 P s lSoP = -9 P s\s lP = -2-8 ± 0.2 ± 0.7 GeV" 1 , 

9t\h 0P = -29t iHiP = -tjl&np = -0-8 ± 0.05 ± 0.2, 

2 



9t\h 0P = 2g d T \ HlP = -2^ -g d T \ Hlp = -9.9 ± 1.3 ± 2.5 GeV" 2 , 
9tIh iP = y| 9t 2 h 0P = ^9t 2 2 h iP = -3.1 ± 0.1 ± 0.8 GeV" 2 , 
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f%s p = 2 9 P t[ SiP = tjl&sv = 35 ± °- 3 ± °- 9 GeV_1 > 
9 P t iSiP = -yjl&sop = -^9%s lP = -3-7 ± 0-3 ± 0.9 GeV -1 , 

9% SlP = &s p = -^9% SlP = 4-1 ± 0.5 ± 1.0 GeV" 3 . (34) 

The second error comes from the uncertainty of /'s. The above relations between coupling constants are consistent 
with the HQET leading order expectation. 

Replacing the p meson parameters by those for the u meson, one obtains the u> meson couplings with the heavy 
mesons: 

9 P h h 0U = ~~9 P H lHlU = -0-29 ± 0.04 GeV 2 , 

~9 P h iHo » = ~9 P H lHl . = -0-41 ± 0.01 GeV 2 , 

~9s h iU = -9s\h 0U = -9s\h iU = -0.36 ± 0.03 GeV 3 , 

9%h xU = -~9s lHoU = -9s\h iU = -0.33 ± 0.05 GeV, 

~9 P s° oSo . = -~9 P s lSl » = -0-29 ± 0.03 GeV 2 , 

9&So U = -9&s lu = -0- 38 ±0 -02 GeV 2 , 

~9t\h . = -2ft lHlU = "2^U 1M = -0.04 ± 0.002 GeV 4 , 
~9t iHoU = *9%h iU = - 2 \H~9tIh iU = -0-43 ± 0.05 GeV 2 , 
9%h iU = \H~9tIh u = yfc9%H lU = -0.13 ±0.01 GeV 2 , 
~9t\s . = 2~9t iSi . = 9% SlU = 0- 15 ± 0- 01 GeV3 - 

9% Slu = -yjl&so*, = -^5 1W = -o- 16 ± o- 01 Gev3 > 

9%s lU = -\H~9t 2 s u, = -^9% SlU = 0.17 ± 0.02 GeV, (35) 



fftoHou, = -9h iHiU = -4.6 ± 0.6 ± 1.2 GeV" 1 , 
9% lHoU = 9% lHlU = -6-6 ± 0.2 ± 1.6 GeV- 1 , 
9s\h iU = -9s\h 0U = -9s\h iU = 3 -6 ± 0-3 ± 0.9, 
9s\h iU = -9i\ HoU = ~9s\h iU = 33 ± 0-5 ± 0.8 GeV" 2 , 
9 P s° s . = -9%s lU = -1-8 ±0.2 ±0.5 GeV" 1 , 
i&Sou, = -9 P s lSl « = -2-4 ± 0.1 ± 0.6 GeV- 1 , 

9t iHo u = -*9%h iU = -tjl&nu, = -0.8 ±0.04 ±0.2, 
9t iHou = Z&Hw = - 2 Vf 9t 3 h iU = -9-0 ± 1-2 ± 2.2 GeV" 2 , 
9% Hlu = \[\9kH 0U = ^9tIh 1u = -2-8 ± 0.1 ± 0.7 GeV" 2 , 
9 P t\s » = *9%s lU = tJl&SiU, = 3 - 2 ± 0-3 ± 0.8 GeV" 1 , 
9%s lU = -J\a^s 0U = -^9% SlU = ~ 33 ± 0- 3 ± 0.8 GeV- 1 
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9t iSi u> = -sjlg f T 2 s . = -^9t 2Si u = 3-6 ± 0.4 ± 0.9 GeV^ 3 . (36) 



IV. CONCLUSION 

We have calculated the light vector meson couplings with heavy mesons in the leading order of HQET within 
the framework of LCQSR. The sum rules are stable with the variations of the Borel parameter and the continuum 
threshold. Some possible sources of the errors in our calculation include the inherent inaccuracy of LCQSR: the 
omission of the higher order terms in OPE, the choice of u> c , the variation of the coupling constant with the Borel 
parameter T in the working interval and the approximation in the light-cone distribution amplitudes of the p meson. 
The uncertainty in /'s and A's also leads to errors. 

The extracted vector meson heavy meson coupling constants may be helpful in the study of the interaction between 
two B(D) mesons. They may play an important role in the formation of these possible molecular candidates composed 
of two B(D) mesons. They may also play a role in the interpretation of the proximity of X(3872), Y(4260) and Z(4430) 
to the threshold of two charmed mesons through the couple-channel mechanism. 
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APPENDIX A: THE p DECAY AMPLITUDES OF HEAVY MESONS 

The definitions of the p coupling constants not presented in the text are 

M(H ^H +p) = (e* ■ q t )g p H ° oHop , (Al) 
M(H, -^H Q + p) = e^^g p H \ Hap , (A2) 
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M{H X ^Ht+p) 
M(S ^So + p) 



M(Si 
M(Si 

M(S - 

M(Si - 

M(Sx - 

M{T X - 

M{Ti - 

M{T X 
M{Ti 



M(T 2 
M(T 2 



50 + p) 

51 + p) 



M(T 2 -> Si +p) 



( e * • 9t)gf oSoP , 

'op' 



(e*-»ft)(e*.g t )-(e*. <)(»,• (&) 



e r,e'gv pl 



► S + p) 

► 5i + p) = (e* • q t )(e* ■ Vt) 9 f lSlP + [(e* ■ ■ ft) - (e* ■ ej i(„ ^ g, I 
ff! + p) = (e* • e t *MU P + [(c* • ft)(e* • ft) - i(e* • 4)ft 2 
H + p) = (e* ■ th)9s\h p + \(v • ft)(e* • ft) - J(e* • ft)ft 2 
Hi+P) 
H + p) 
Hi+p) 



pi 

9h 1 h 1p > 



n dl 
n dl 



e Vt'<lv( e *. qt} __ e m*e'v q 2 



9s 1 H 1 p> 



(e* • m)9T lHoP + [(V ■ ft)(e* ' ft) ~ g(e* • Vtht 



9t 1 h oP -> 



^ V 9 6 t\h 



IP 



di 
9t 1 h iP 



J (ft-0 + e £ eqv (qt-v) 



di 

9TiHip> 



■qe*qv pi 
6 9t 1 S oP > 



(e* • g t ) -(e* •£?)(»/. ft ) 



pi 
9t 1 s iP 



(e* t ■ »?)(ft ' O + (ft ' ^)(e t * • £*) - g(e t * • r?)(e* • g t ) 



p2 

9t 1 s iP 



+ {(ft-^)(ft-e*)(e*-ft) 

'fat • e*)(e* ■ q t ) + (e* t ■ V )(q t ■ e*) + (q t ■ V )(e* t ■ e*)] }g% 



Q[_ 

5 



Sip' 



#o + p) 

H l+ P) = Va ia2 



r}a ia2 \ e Qt + 6 ft )9t 2 h qP i 

2 . 



- 3 9t 



(e*t ■ e*) 



9t 2 h iP 



~\~Vaia 2 
1 



{ [e* Ql 9? 2 + e* Q2 ft Ql - g5r ia2 (ft ' O] ( e * • ft) 



\9tIh iP 



+ry QlQ2 {2[er i ft a2 (ft • O + ft ai G r 2 (ft ■ O " 2ft Ql ft tt2 (e t * • e*) 
e *a lq « 2 + e *«* q ai _ 2 g^( qt . e *)] (e* • ft) 

2 9 r a2 (e* t -e*)]q?}g% HlP , 



M(T 2 -> 5 + p) = 7y aiQ2 



! (e*-ft) 



P2 

9t 2 s oP 



+ Vaia2 {q^qr(e* ■ q t ) - f [^(e* • ft) + e^g* + g^ 2 ] } 



/2 

9t 2 s oP ' 



1Ja 1 ct 2 
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9t 2 Si P 
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ft 
5 



,aie q-u *«2 _|_ o 2 e qv *ai . ai£ e u o 2 . a 2 e e u ai 



'7/ 



ft 



p2 

9t 2 s iP 



(9t 2 Si P - 



(A3) 
(A4) 
(A5) 
(A6) 

(A7) 

(A8) 

(A9) 

(A10) 

(All) 
(A12) 



(A13) 
(A14) 



(A15) 



(A16) 



(A17) 



Note that these decay amplitudes may be organized in another way. For example, the tensor structure corresponding 
to g P H 1 H lP was defined as (e* • v)(e* ■ -qt) in equation (28) of Ref. § rather than (e* • ft)(e* • f] t ) in Eq. (|A3[) . Since we 
have (e* • ft) = —(q ■ v)(e* ■ v), the essentially same sum rule as Eq. (24) of Ref. Q can be obtained if we isolate the 
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tensor structure (e* • v)(e* -ft). 

To derive sum rules for these coupling constants, 



/ 



/ 



d 4 xe- ikx (p(q)\T{J 0t _ th (0)jl ,(x 
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}|0> 
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}|0) 



we consider the following correlators: 
+ (erq?~q?ef)G P l Hip (co,to>), 



(A18) 
(A19) 



(A20) 
(A21) 

(A22) 
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}|0> 
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di 
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-( ff f l/3 er 2 +5 t Q2/3 er i " 2 9 r a2 ef )q 2 t 
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,a 1 0e'v^a 2 _|_ £ a 2 0e'v ^ai _|_ ^aifiqv e *a 2 _|_ ^a 2 fiqv g*«l 

L ai ^ v qf 2 (e* ■ q t ) + t a2l3qv q^{e* ■ q t ) 
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^a-iPqv *a 2 , a 2 flqv *a t , a 1 fle*v a 2 , a 2 f3e*v ol\ 



Qt 



Qt 
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APPENDIX B: THE p MESON LIGHT-CONE DISTRIBUTION AMPLITUDES 

The definitions of the distribution amplitudes used in the text read as [l3[ 



(0\u(z)j p d(~z)\p-(P,\)) = f p m p 



P P — \ due^ z ^(u^ 2 ) +e^l / due^ z g { l\u^ 2 ) 



P-z Jo 
1 e W ■ z 

2 Zfi (p-zr-""j 



m 2 / due^ z g 3 (u,n 2 ) 



(0\u(z) W d(-z)\p-(P,X)) = \f pmp e^e { H Pa zp f due^g^ (u, 



(Bl) 
(B2) 



(0\u(z)<j p „d(~z)\p-(P,\)) = if p 



+ {p p z v -p v Zfi)j r^TOp / due lip ' z h^(u,n 2 ) 

' 1 ' o 

,2 ,.1 



(p ■ zf 



(0\u(z)d(-z)\p-(P,X)) = -tfJ(eWz)m 2 p / du e^ p ' z h^\u, p 2 ) 

Jo 



The vector and tensor decay constants f p and fj are defined as 



|u(0) 7 ^(0)|p-(P,A)> = f p m p el x \ 
fi(0)^d(0)|p-(P,A)) = z/J^P.-e^). 



The distribution amplitude 0|| and <^>j_ are of twist-2, g^\ 9±, foji and ^z. ^ & rc twist-3 and 33, /13 are twist-4. All 
functions <f> = 4>±, g^\ g^\ g 3 , h^} are normalized to satisfy Qdu<j){u) = 1. 

The 3-particle distribution amplitudes are defined as [l3| 

(0\u(z)gG lll/ j a j 5 d(-z)\p'(P,X)) = f p m p p a {p„e'i ) - p p e ( ll]A(v,pz) 
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+fj™? p \p a e { llgj v - Ppe^g^ - p a e { Hg^ + Ppe^g^T^ (v,pz) 
+fjm 2 p \p»e ( £l9l3v - P^%9a V - P^ ) a gp fl + 2Ve^s^jT 2 (4) \v,pz) 
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pz 

Til? 
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pz 

+ ... (B9) 

(0\u(z)gG^(vz)d(-z)\p-(P,X)) = i.fjm 2 p {e { iy - e { Hp,}S(v,pz), 
(0\u(z)igG^(vzh 5 d(-z)\p-(P,X)) = if p ml\e%p v - e™ Pli ]S(v,pz). (BIO) 



where 




(Bll) 



etc. The integration measure is 




(B12) 



The distribution amplitudes A, V and T are twist-3 and the others are twist-4. 



